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We develop the perturbative QCD formalism for inclusive heavy hadron decays. 
Transverse degrees of freedom of partons are introduced to facilitate the factoriza- 
tion of the heavy hadron decays. 



In this talk, we shall derive the PQCD factorization formula for the semilep- 
tonic decay B X„£i/ in a rigorous way, to demonstrate how to use techniques 
of PQCD for heavy hadron decays. The factorization procedures demand the 
inclusion of the transverse degrees of freedom of partons at the end point of the 
charge lepton spectrum where an outgoing jet is present. Hence, we have to 
perform the resummation of large perturbative correictions in the transversal 
configuration space using the technique developed inB, which is also accurate 
up to next-to-leading logarithms. 

We work in the rest frame of the B meson, and choose the following light- 
cone components for relevant momenta, 

Pb = {P^,Pb,0±) , P£ = (p,+ ,0,Oi), p. = (p^,p;,p.±) , (1) 

with = Pg = Mb I and = 0. The independent variables are identified 
as p/, p~ and pj, and their relations to Eg, and qo are Ei = p^/-\/2, 
q^ = 2p^p^ and = (p/ +pt +P^)/V^, respectively. We define Pb = Pb -p 
as the b quark momentum, which satisfies P^ « M^, Mj, being the b quark 
mass, p is the kicks from the light components inside the B meson, which has 
a large plus component p"*" and small transverse components p_L . The b quark 
decays into a u quark with momentum P„ = Pb —p — q. We have distinguished 
the B meson momentum Pb from the b quark momentum Pi, here. 
It is more convenient to employ the scaling variables 

■2Ee (?2 2qo 

instead of the dimensionful ones En, q^ and q^. Note that the scaling variables 
are defined in terms of the B meson mass Mb, since we formulate the factor- 
ization according to the B meson kinematics. For massless leptons, it is easy 
to show, using the momentum configurations defined in eq. (||) , that the phase 
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space is given by 

0<a;<l, 0<y<a;, -+x<yo<y + l. (3) 

X 

In the end point region with x ^ 1 {p^ and y — > {p~ — > 0), 

we have j/o ~^ 1 {pt ^ 0) ^i^^ p ^ 0. The u quark then has a large minus 
component = (1 — y/x)MB/^ but a very small plus component F+ = 
(1 — 2/0 — ylx)MB/V^^ and thus a very small invariant = — yo + y)j 

which forms an on-shell jet subprocess. The u quark travels a long distance 
of ©(I/Aqcd) before hadronized. Besides, the B meson is dominated by soft 
dynamics, which is the origin of the soft function stated in the Introduction. 
The remaining dominant subprocess is the hard one, which contains the weak 
decay vertices. Therefore, the important contributions are factorized into the 
soft (5), jet (J) and hard {H) subprocesses. 

The factorization formula for the inclusive semilcptonic decay B — > X^iv 
is written as 

xS{z,p±, fi)j{z, P.;; ,p±, fi)H{z, ,p±, n) , (4) 

with the momentum fraction z defined by z = P^ / Pb — 1 ^ / Pb ^'^'^ 
r^""* = j^^lKibpA^B- M sq. (Q) is the renormalization and factorization 
scale. The triple differential decay rate is, of course, /j, independent. Note that 
in the region y ^ x ^ 1 the outgoing u quark becomes soft and eq. (Q) fails. 
We shall show that contributions from this dangerous region are suppressed by 
phase space. The upper limit of z takes the value Zmax = 1 in our analysis. If 
performing the factorization according to the h quark kinematics, one must as- 
sume Zniax = Mb /Mb, which is greater than 1, in order to fill up the kinematic 
window. It has been explained □ that Zmax > 1 is not allowed in perturbation 
theory, and is thus of nonperturbative origin. From the kinematic constraints 
in eq. and the on-shell condition of the w-quark jet, the lower limit of z 
should be Zmin = X, instead of z,„i„ = 0. 

The tree-level expressions for the convolution factors J and H are given 

by 

J(°) = S{P^)^S^Ml 

cx {Pb-P.){pi-Pu) 

oc {x-y){yQ-X' 



l-yo + y~{l-z){l~^) 

X 

= {{Pb -p)-p.){PfPu) 
^ X Mi 



2p± ■ 

Ml 



Ml 



(5) 
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Equation can be regarded as an expression at the intermediate stage in the 
derivation of conventional factorization theorems. If the pj_ dependence in J 
and H is neghgible, the variable in S can be integrated over, and eq. (||) 
reduces to the conventional factorization formula. However, it is obvious from 
eq. (H) that at least the pj^ dependence in J is not negligible, especially in 
the end-point region. This is the reason we introduce the transverse degrees 
of freedom into our analysis. 

Suppose we consider higher-order corrections to eq. (Q) from a gluon cross- 
ing the final state cut, and route the loop momentum t through, say, the jet sub- 
process. Without losing generality, we approximate J{p^ +^^, Pu +^ i P± + 
t^) K, J{p'^,P~,p± + £±) according to the kinimatic relations £^ < p'^, 

< and £± « pj_. Hence, the loop integral cannot be performed un- 
less the dependence of J on transverse momentum is known. This difficulty 
can be removed by Fourier transform, 

J(p+,p-,p^+^^) = J ^jV,p-,b)e^(P-+^-) \ (6) 

where the impact parameter b (Fourier conjugate variable of p±) measures the 
transverse distance travelled by the jet. Using eq. (||), the £± dependence is 
decoupled from the jet function, and the factor e'^^-^'^ is absorbed into the loop 
integral, which can then be performed. Therefore, an extra factor e^^^'^ is 
associated with each gluon crossing the final state cut in our formalism. 

To further simplify the factorization formula, we neglect those terms in- 
volving p^x in J and H. This is a good approximation for a; — > and x — > 1, 
since for a; — > contributions from transverse momenta are not important, 
and for a; ^ 1 the magnitude pi,± = M b \/ {yo — x — y/x)y/x vanishes. Equa- 
tion (H) then becomes 

1 d^r _ 2 f\ f 
^d^d^o ~ ^^^y. J (2^ 

x^(z,b,^)J(z,P„",b,/i)i?(z,P„~,^) . (7) 

It can be shown that the dominant subprocesses J, S and H contain large 
logarithms form radiative corrections. In particular, J gives rise to double 
(leading) logarithms in the end point region. These large corrections spoil the 
perturbation theory and must be organized. One has to resum these large 
corrections up to next-to-leading logarithms. The first step in resummation is 
to map out the leading regions of radiative corrections. In the coUinear region 
with the loop momentum £ parallel to Pu and in the soft region with £ ^ we 
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can eikonalizc the heavy b-quark hne. With the eikonal approximation, the fa- 
quark propagator is expressed as 1/v ■ i to the order 1/M b with w = (1, 1, Oj^). 
Hence, the factorization in eq. (0) is in fact vaUd to 0{\IMb)- The physics 
involved in this approximation is that a soft gluon or a gluon moving parallel 
to Pu can not explore the details of the b quark, and its dynamics can be 
factorized. This is consistent with the HQEFT, where the h quark is treated 
as a classical relativistic particle carrying color source. Since the large mass 
Mb does not appear in the eikonal propragator, the only large scale in J is 
P~ . After eikonalization procedures, one obtain the Sudakav exponent that 
resumes large double logarithms arise from both coUinear and soft regions. 



Jl/h p Jl/b 



A = C, 



67 
~9 



TT 

y 



10 



27 



(8) 
(9) 



with Cf = 4/3 the color factor. 

Having summed up the double logarithms, we concentrate on the single 
logarithms in S, H and the initial condition J{b, fi). Since both the differential 
decay rate and the Sudakov exponent s{P~,b) are RG invariant, we have the 
following RG equations: 



with 



VJ{b,^l) = -27, J(6,/i), 
VS{b,^) = -'^sS{b,fi), 
VH{p-,^i) - (27, + 7s)i/(F-,M) 



(10) 



(11) 



7q = — tts/TT is the quark anomalous dimension in axial gauge, and 75 = 
— (as/7r)CF is the anomalous dimension of S. Integrating eq. (|lO|), we obtain 
the evolution of all the convolution factors, 



Jiz,p-,b,fi) 
S(z,b,fi) 
H{z,p-,ii) 



exp 



exp 



exp 



-2s(p-,6)-2 r ^7,(«.(m)) 



J>,5,l/6), 



^7s(as(Ai)) 
1/6 M 



f{z,b,l/b), 



■27,(a,(/2))+7s(a,(/2))] 



Hiz,p-,pm 
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We shall neglect the intrinsic b dependence of the distribution function / be- 
low. If Sudakov suppression in the large-6 region is strong, we may drop the 
evolution of / and J in 5, which is proportional to as{l/b). Hence, we as- 
sume f{z,b, 1/b) = f{z), J{z,b, 1/b) — J^^\z,b), the Fourier transform of the 
tree-level expression in eq. (||), and H{z, P~ , P~) = H^^\z,P~). 

Substituting eq. into (0), we derive the factorization formula for the 
inclusive semileptonic B meson decay, 

' '^'^ - Mil'dz r'^nz)ji^H^m'"H^,Pu-) 



rfUxdydyo "J, J a 27r ' 

xexph^(p-,6)] . (13) 

The complete Sudakov exponent S is given by 

S{P-,b) ^ 2s{p-,b) ~ T^ln^ , (14) 
3pi b 

with P^ = ln(Pj7/A), which combines all the exponents in eq. (|l2|) and in- 
cludes both leading and next-to-leading logarithms. It is straightforward to 
observe from eq. (|lj) that the Sudakov form factor falls off quickly at 
large b ~ 1/A, where as{l/b) > 1 and perturbation theory fails. Hence, the 
Sudakov form factor guarantees that main contributions to the factorization 
formula come from the small b, or short-distance, region, and the perturbative 
treatment is indeed self-consistent. We stress that our formalism is applicable 
to the entire range of the spectrum. 

Since all the double logarithmic corrections have been absorbed into the 
jet subprocess, the soft function S contains only soft single logarithms. These 
single logarithms can be summed by solving the RG equation VS = —jsS as 
shown in eq. (10). The solution has been given in eq. (p^), 



S{z, 5, fj.) — exp 



^7s(as(Ai)) 
1/b M 



f{z,b), (15) 



where the initial condition f{z) for the evolution of S must be determined 
phenomenologically. From the definition of S, it is obvious that / depends 
only on the properties of the bound state \B), but not on the particular short 
distance subprocess. Therefore, / is a process-independent universal function 
describing the distribution of the b quark inside a B meson. Since the intrin- 
sic b dependence (the perturbative b dependence has been collected into the 
exponent S) is not known yet, we take the ansatz f{z,b) = /(z) exp(— E(6)), 
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which leads to S ^ as 6 according to the definition /(z, 6 = 0) = f{z). 
It is also natrual to assume S > for all b from the viewpoint that the b 
quark is bounded inside the B meson. Hence, the intrinsic b depedence pro- 
vides further suppression. The nonperturbative function /(z), identified as the 
_B-meson distribution function, can be expressed as the matrix element of the 
6 quark fields, whose first three moments areoij 



f{z)dz = 1 , 







f{z){l-z)dz^ — 
niB 

mi 



+ 0(A|cdMI) 

Ai 
3 



/(z)(l - zfdz = _ ( - ^ ) + OiAl^cu/ml) 



(16) 



Though the exponent E is unknown, we can, however, extract its leading 
behavior by means of the IR renormalon analysis. Note that the perturbative 
Sudakov factor in Eq. (|l^) becomes unreliable as the transverse distance b 
approaches 1/Aqcd- Near this end point, as{l/b) diverges, and IR renormalon 
contributions are significant. We reexpress the RG result of the evolution of 
the distribution function, which is contained in the second term of S, as 



W — exp 



47rC, 



(27r)4 (v ■ I) 



^^^^27r^(Z2)a,(;2)e*i-bArM^ 



(17) 



The loop integral corresponds to the correction from a real soft gluon attaching 
the two valence b quarks, whose propagators have been replaced by the eikonal 
lines in the direction = (1, 1,0) Q. The tensor Nt'" = g'"' - {n'^l'' + l^'n'')/ {n- 
I) + n?l^l'^ / (n • /)^ comes from the gluon propagator in axial gauge n ■ A — 0. 
We have set the argument of the running ag to l"^, which is conjugate to the 
scale b of the distribution function. 

Substituting the identity as{l'^) = tt Jg rfcr exp[— 2(t/3i ln(/T/AQCD)], Pi — 
£33 — 2n/)/12, into Eq. (p7[), and performing the loop integral for n cx (—1, 1, 0) 
B, we obtain 



W — exp 



da 



^bA. 



QCD 



2cr/3i 



r(-a/3i) 
r(i + a/3i) 



(18) 



It is easy to observe that the pole of r{—a(3i) at ct ^ gives the perturbative 
anomalous dimension of the distribution function appearing in Eq. (^). The ex- 
tra poles at a l//3i, 2/ Pi,..., then correspond to the IR renormalons, giving 
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corrections of powers 6^, 6^,..., respectively. These renormalons generate singu- 
larities, which must be compensated by the nonperturbative power correctons 
in order to have a well-defined perturbative expansion. Using the "minimal" 
ansatzH, ie picking up only the leading (left-most) renornialon contribution, 
we parametrize the exponent S(6) by S(6) = c'm^b^, corresponding to the 
fact that the power corrections start at 0{lP'). Certainly, other parametriza- 
tions consistent with this fact are equally good. We postulate the following 
two-parameter distribution function for the B meson, 

fB{z) = N 6{l - z) , (19) 

[{z — a) + ez\^ 

The charged lepton spectrum for the decay B — )■ Xuf-v from the naive 
quark model is obtained by simply choosing f{z) — 5{1 — z) and ignoring the 
transverse momentum dependence in j'"' and H'^^\ A simple calculation leads 
to 

1 

(3 -2a;) . (20) 



r/o) dx 6 

This spectrum does not fall off at the end point of the spectrum, contradicting 
the observed behavior of the inclusive semileptonic decays of B mesons. The 
discrepancy implies that the tree-level analysis is not appropriate, especially 
in the end-point region where PQCD corrections are important. 

The spectrum from the parton model without Sudakov suppression is ob- 
tained by adopting H^^^ — (x — y)\yQ — x — (\ — z)y/x] and = Mg[l — yo + 
y — {1 — z){l — y/x)]. With integration over j/g, we derive 

1 dr /"^ 

dy / dzf{z){x^y){y + z-x) . (21) 



r/'') dx 



At last, including Sudakov suppression and all other single logarithms, we 
arrive at the charged lepton spectrum of the B X^tv decay that takes into 
account both large perturbative and nonperturbative corrections. 



1 dr r f^/^ 

- Mb dy db dzf{z){x^y)^ 

Jo Jo Jx 

-^J2{^MBb)+^^M^MBb) 



r/o) dx 



MBb 



{z + y-x)Ji{^MBb) 
siP-.b) ^ (22) 



with ^ = {x — y){z/x — 1) and Ji, J2, ^3 are the Bessel functions of order 1,2 
and 3 respectively. 
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We stress that our results do not violate the conclusion from HQEFT, 
if they were interpreted in a proper way. To confirm this, we identify Pf, — 
{M^ /2Pg , Pg , 0) as the momentum carried by a free b quark in factorization 
theorems for B meson decays, where the minus component Pf^ has been set 
to Pg . That is, the free b quark is not at rest inside the B meson. We then 
reexpress eq. ( ^T| ) into a form similar to that ind: 

I dT , , / \ fM^\ , , 



with F{x) = x^(3 — 2x)/6 being the quark-model prediction derived from the 
conventional approaches and 



F[jit)M{x)^ X "^yj dzf{z)ix-y)iy + z-x)-F{x)9i^j^ 

/ _ ^ . . ^^^^ 

The step function in eq. (23) specifies the maximal Eg in the decay of a free b 

quark with the above momentum Ph. The function M{x), representing nonper- 

turbative corrections to the b quark decay, coincides with the shape function 

S{x) defined incl. 

We shall show that the contribution from M (x) to the total decay rate is 
indeed of 0(1/A/^). Integrating eq. ( |2^ ) over x, we obtain 

j^j M{x)dx = - dzz^fiz) - F{x)dx . (25) 

An arbitrary structure function /, which possesses the same moment as in 
eq. (p^, can be expanded in terms of (5- functions: 

f{z) = 6(1 -z)- ^5' {I -z) + 0{K^/Ml) . (26) 

Inserting eq. ( p6[ ) into (^), we justify straightforwardly that the nonperturba- 
tive correction 

Mb\ ^ , , , 1 A 



FItT ^^(^)^^ = / ^(^)^^-T^7^/ dzz^S'il~z) + OiA'/M's) 

\MbJ J Ml /Ml 12 Mb J Q 

vanishes at 0{1/Mb) as concluded ini. In summary, the quark-model con- 
tribution from the window between x = M^/M]^ and x = 1, with a width of 



0{1/Mb), cancels the 0{1/Mb) correction from the structure function, such 
that the nonperturbativc correction is of ©(l/A'/g). Our result provides an 
explicit dynamical demonstration the validity of global quark-hadron duality. 
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